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ABSTRACT 

We have discovered a dwarf nova (DN) of type SU UMa in Kepler data which is 7.0 arcsec 
from the G-type exoplanet survey target KIC 4378554. The DN appears as a background 
source in the pixel aperture of the foreground G star. We extracted only the pixels where 
the DN is present and observed the source to undergo five outbursts - one a superoutburst 
- over a timespan of 22 months. The superoutburst was triggered by a normal outburst, a 
feature that has been seen in all DNe superoutburst observed by Kepler. Superhumps during 
the superoutburst had a period of 1. 842 ±0. 004 h and we see a transition from disc-dominated 
superhump signal to a mix of disc and accretion stream impact. Predictions of the number of 
DNe present in Kepler data based on previously published space densities vary from 0.3 to 
258. An investigation of the background pixels targets would lead to firmer constraints on the 
space density of DN. 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are binary star systems consisting 
of a white dwarf primary and a less massive main sequence or 
post-main sequence secondary. The two stars orbit a common cen- 
tre of mass with typical periods from between 50 min and 5.7 d 
jDownes et alj|200ll ; iRitter & Koibll2003l : iThorstensen et alj|200fj| : 
Seki guchilll992TT The small physical size of the binary system 
causes the secondary to be distorted by the gravitational field of 
the white dwarf and leads to the filling of its Roche lobe. Mass 
is transferred through the LI inner Lagrangian point and forms a 
stream of material falling into the potential well of the white dwarf. 
In the case where the white dw arf is non-magnetic, the inf a iling ma- 
terial forms an accretion disc Jshakura & S unvaevlll973h . Viscous 
shearing layers amplify a weak magnetic field within the disc, and 
lead to turbulent angu lar momentum transport radially outwards 
tealbus&Hawldll998h . Mass can then move inwards and be ac- 
creted onto the surface of the white dwarf ( Helliei 2001). At optical 
wavelengths, the disc is usually much more luminous than the white 
dwarf, its companion, and the impact of the accretion stream onto 
the disc. 

The DN su bclass of CVs undergo quasi-periodic outbursts 
(see|Warner 2003, for a review). During an outburst the optical 
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bright ness of the system can increase by a factor of 10 or more 
(e.g. lRamsavetalj2009l) . which interrupts the faint quiescent state 
in which the DNe typically spend most of their time. 

If the mass transfer rate from the secondary to the accretion 
disc is greater than the rate of disc accretion onto the white dwarf, 
M(2) > M(d), then the disc radius will increase over time. This 
can only happen if the dimensionless disc viscosity parameter a 
is low (where low is usuall y taken in models to be < 0.03, e.g. 
ISmaklll984 ICannizzo|[T993n . The disk instability model predicts 
that eventually the surface density in the disc will exceeds a crit- 
ical v alue S max , which initiates a phase transition in the disc ma- 
terial ( Hosh| 1 19791 ; lLasotal 12001^ At this point the a param eter 
increases (by a factor of ~ 10 ISmaklll984l ; ICannizzolll993h and 
M(2) < M(d) for a time until the surface density at the outer edge 
of the disc decreases below some critical surface density E m i„, and 
q returns to the original lower valu e. This model is known as the 
limit cycle accretion disc instability (Lasota 200l]). 

A subset of DNe known as SU UMa stars occasionally un- 
dergo a superoutburst (SO) in addition to undergoing normal out- 
bursts (NO). SOs are ty pically a facto r of two brighter and 5-10 
times longer than a NO jWarnerlliooH) . Superhumps, photometric 
variability on a period a few per cent longer than the orbital pe- 
riod of the system, are observed du ring the SO. During a NO the 
disc expands in radius (Osaki 1989). After a number of these out- 
bursts, the disc can become large enough to grow beyond the inner 
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Lindblad resonance, close to the 3:1 co -rotation resonance with the 
binary orbital period! Whitehurst 1988). Once there is enough mass 
at this resonance radius, the superhump oscillation can be driven by 
tides. 

We identified a background source in Kepler spacecraft 
archival data whose light curve resembled that of a DN. In Sj2]we 
describe the observations made by the Kepler spacecraft of this 
source, our re-extraction of the light curve from pixel level data 
using a more appropriate pixel aperture, and our removal of sys- 
tematic photometric effects from the flux time series. In |j3]we dis- 
cuss observed outbursts, and in Sj4]we discuss the implications of 
these observations on the nature of outbursts from DN. 



2 OBSERVATIONS 

Kepler has been almost c ontinuously observing over 1.6 x . 10 5 st ars 
at a cadence of 30 min teorucki et alj|20ict iKoch et alJkoid) at 
milli- to micro-magnitude precision IGilliland et alj201 If) since the 
spacecraft's launch on 2009 March 06. While the primary aim of 
the Kepler mission is the discovery of exoplanets, the long base- 
line of continual, high precision observations has also provided an 
unrivalled data-set to the stellar astrophysics community. 

The available bandwidth does not allow all starlight falling 
onto the Kepler detectors to be downloaded at a cadence of 30 
min. Instead the pixels to be collected must be preselected and 
account for only 5.75 per cent of the total detector area available 
fcoch et alj|2010h . 

The continuity of observations are broken every month by the 
satellite pointing towards Earth for data to be downloaded. Once 
every three months the satellite rolls 90° to keep the sunshade 
pointed in the appropriate direction but has the effect of changing 
the CCD that a star falls upon every quarter. Because of the change 
in physical position of the instrument and because a new data-set is 
downloaded, the data from Kepler is divided into quarters starting 
with Quarter (Q0), containing commissioning observations, and 
full science operations beginning in Ql. 

Each observed star is extracted by selecting a pixel mask 
which can change in size and shape every quarter because the po- 
sition of the star in pixel space is different. For this reason the con- 
tamination from background stars is different every quarter - this 
contamination can be significant because part of the Kepler field 
covers the Galactic plane (where the stellar density is high) and the 
pixels are relatively large (3.98" x 3.98"). 

The observations described in this paper are of one particular 
target aperture - that around KIC 4378554. This star was chosen 
to be a target in the Kepler exoplanetary program. Kepler observed 
KIC 4378554 during Ql-8 excluding Q6 (see Table|T]for a log of 
the observations). No planet has been found transiting this system 
as yet. The Kepler Input Catalogue lists KIC 4378554 as a K p = 
14.9 sta£|with colour (g'-r') = 0.57, an effective temperature of 
5617 K, and a surface gravity of log g — 4.615 in cgs units; these 
valu es are consistent wi th the source being a main sequence G-type 
star jBrown et alj201 lb . The star is only visible during three of the 
four seasons of observations because it lies on Module 3 for one 
quarter. Module 3 failed during Q4. 



1 Kp refers to approxi mate Kepler instrum ental magnitudes. For more de- 
tails of this passband see Brown et al. (201 1). 



2.1 Optimal aperture observations 

Two light curve products are provided by the Kepler project: light 
curves extracted using Simple Aperture Photometry (SAP) and 
light curves which have underg one the additional s tep of Pre- 
Search Data Conditioning (PDC. iTwicken et al.ll201oh . The latter 
of these is produced for use in exoplanet searches but is known to 
cause problems if used for stellar astrophysics work. These prob- 
lems include the attenuation of astrophysical signal and the injec- 
tion of spurious signals, because of this we use the SAP data. A 
light curve showing Q5 data can be seen in FigQ] These data have 
not had systematic artefacts removed, and there are several notice- 
able features due to systematics such as the increase in flux level 
over the final 50 d of the quarter and the thermal settling after Earth- 
points at BJD 2 455 307 and 2 455 337. However, there is one ob- 
vious feature that appears to be astrophysical in nature - a 2 per 
cent increase in brightness with a peak at BJD 2 455 342.8 lasting 
2.5 d. In addition to the brightening shown here, we identified four 
similar events - two in Q3, one in Q4 and one in Q7. Events like 
these are not typically seen in G-type stars and therefore warranted 
further investigation. 

2.2 Custom light curve extraction from target pixel files 

In order to determine whether the brightness increases are due to 
KIC 4378554 or a background star in the optimal aperture of KIC 
4378554, we examined the archived target pixel filefl These files 
contain a series of images - one for each time-stamp - showing the 
individual pixels downloaded from the Kepler spacecraft. These 
images were calibra ted by the Kepler pipeline (as described by 
Ouint ana etafl feoiO). Shown in Fig. [2] is the Q5 time-series for 
each calibrated pixel in the aperture of KIC 4378554. The pixels 
shown in grey were used to create the light curve in FigQ] From 
Fig.[2]we can see that the brightness increase at BJD 2 455 342.8 
is not caused by KIC 4378554 but from a source which appears to 
have its centre on or close to pixel (1043,857) approximately 1.5 
pixels south-west of the target star. For convenience we will refer 
to this source as NIK 1 (standing for Not In KIC). 

We extracted a new light curve by summing only those pix- 
els which NIK 1 contributes flux using the KEPMASK and KEPEX- 
TRACT tools from the PyKE Kepler community softwar^fl In quar- 
ters where we see a sudden brightening event, this is relatively 
straight forward. This is the case in Q3, 4, 5 and 7. For Ql we use 
the Q5 aperture and for Q8 we used the Q4 aperture because NIK 1 
will be at the same position on the CCD four quarters apart. For Q2 
we estimate the position of NIK 1 from the full-frame image taken 
at the end of Q3, when a brightening event occurred and NIK 1 is 
visibl^B We overlayed the full-frame image containing NIK 1 on 
the pixel level data and selected the appropriate pixels. The number 
of pixels we use for all quarters is listed in Table|2] 

2.3 Removing systematic trends using cotrending basis 
vectors 

The light curves we extracted from target pixel files still contain 
systematic features. This is because the extracted aperture does not 

2 Kepler target pixel and light curve files are distributed at 

|http : //archive ■ stsci ■ edu/kepler| 

a PyKE is available from |http : //keplergo . arc ■ nasa ■ gov/ContributedSof twa": 
4 Kepler full-frame images are collected once a month and are available 
from http : //archive . stsci . edu/kepler/f f i/ 
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Table 1. Log of Kepler observations of KIC 4378554 up to BJD 2455635. The start and end Barycentric Julian Date (BJD) refer to the mid-time of the first 
and last cadences on which there are no flags questioning the quality of the data. We find the data quality flags in the light curve and target pixel files under 
the column SAP.QUALITY - non-zero indicates a potentially compromised cadence (the specific meaning of each flag is given in the Kepler Archive Manual 
Fraquelli & Thompson 201 1). Likewise, the number of cadences column in this table excludes cadences which contain potentially compromised data. 



Quarter 


Start BJD 
-2454900 


End BJD 

-2454900 


# cadences 


Module 


Output 


Channel 


1 


64.5115 


97.9832 


1611 


15 


3 


51 


2 


102.7648 


191.4672 


4048 


3 


3 


7 


3 


193.2449 


282.4947 


4099 


11 


3 


35 


4 


285.3758 


375.2017 


4092 


23 
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79 
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376.5095 


471.1621 


4460 


15 
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51 
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Figure 1. The Q5 Simple Aperture Photometry (SAP) light curve made using the optimal aperture calculated for KIC 4378554. The light curve contains 
systematics, such as differential velocity aberration, which causes the change in the general flux level over the last 50 days of the quarter, and a data-gap from 
Earth-point followed by a thermal settling at BJD 2 455 307 and then again at BJD 2 455 337. In addition, there is a sudden, temporary 2 per cent increase in 
the flux level, reaching a peak at BJD 2 455 342.8. 



Table 2. Pixels used to extract the light curve of NIK 1 in the optimal and custom apertures, the time spans masked out of the cotrending basis vector fit and 
the number of basis vectors used in these fits. 



Quarter 


Pixels in optimal ap. 


Pixels in custom ap. 


Time range masked in CBV fit 
BJD - 2454900 


# of CBVs in fit 


1 


6 


2 




3 


2 


5 


1 




3 


3 


6 


5 


272.71465-282.50147 


5 


4 


4 


2 


285.35980-293.98626 


5 


5 


5 


2 


441.89137-448.34535 


6 


7 


6 


5 


621.90676-629.78812 


6 


8 


4 


2 




3 



© 0000 RAS, MNRAS 000, 000-000 



4 Barclay et al. 



< 




i — i — i — i — 1= — i — i — i — i — == — i — i — i — ^ 




1041 



1042 



1043 1044 

Time 



1045 



1046 



Figure 2. Plotted is the flux time series within each pixel collected for KIC 4378554 during Quarter 5, where a Quarter is the number of days between 
each spacecraft roll. The quarter length is typically 91 days but can vary because of, for example, spacecraft anomalies. The numbering around the 
plot provides the pixel column and row on the CCD channel. The blank boxes contain pixels that were not within the pixel mask. The five boxes with the 
grey background are those within the optimal aperture of KIC 4378554 and were used to make the plot shown in Fig[T]and the archived light curve at the 
Multimission Archive at STScI (MAST). The brightening event is not coming from the source within the optimal aperture but from another source on the 
bottom edge of the pixel mask. The flux range is scaled separately for each pixel, the amount of scatter is in each pixel is indicative of the signal-to-noise of 
the data in that pixel. Differential velocity aberration causes KIC 4378554 to move towards the lower left corner over the length of the quarter. The thermal 
settlings cause KIC 4378554 to move from left to right. 



contain all the flux and also contains contaminating flux from the 
target star, KIC 4378554. The amount of flux in the aperture, both 
from NIK 1 and from other contaminating sources, changes contin- 
ually over time due to differential velocity aberra tion and suddenly 
due to changes in the focus after Earth-points l lChristiansen et al] 
1201 lah . 



The instrumental systematics which affect a single light curve 
are unique. However, we see strong correlations between stars 
which are on the same CCD channel. We can take those stars which 
are most highly correlated, which we call the reference ensemble, 
and identify the correlated features. These correlations can be rep- 
resented as linear combinations of orthogonal vectors, which can be 
fit and subtracted from each individual light curve to remove instru- 
mental features. The Kepler team have created a product containing 
these vectors, known as cotrending basis vectors (CBVs), and have 



made them availa ble to the publijf] (a description of how they we re 
created is given in lChristiansen et aU201 lbTlBarclav et alj|201lh . 

The available files consist of the 16 vectors, of length the same 
as the light curve, which provide the greatest contribution to the 
variance of the reference ensemble. There is one set of 16 vectors 
per quarter for every CCD readout. It is crucial that an appropriate 
number of vectors are chosen to fit to to the light curve data; using 
too few will leave instrumental signal in the light curve, using too 
many will cause over-fitting and remove astrophysical signal. 

We used the KEPCOTREND tool from PyKE to perform a lin- 
ear least squares fit of the CBVs to the light curve we extracted 
using KEPEXTRACT, and masked the regions where the brighten- 
ing events occur out of the fit because these regions are dominated 
by astrophysical signal not systematics. We selected the number of 



CBVs are available from http : / / archive . stsci . edu/kepler /cbv . html 
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CBVs to fit by experimenting with the fit: we assumed that there 
was no intrinsic long period variability outside of the brightening 
and attempted to remove any long term trends. Starting with two 
CBVs, we iteratively increased the number of CBVs until the re- 
gions outside of the brightening events was approximately constant. 
The number of CBVs used and the regions masked from the fit are 
listed in Table[2] 

The CBVs only correct for the amplitude of systematic effects 
and not for constant offsets caused by differences between the con- 
tamination from nearby sources and changes in the fraction of flux 
of NIK 1 in the pixel aperture in different quarters. After removing 
the instrumental systematics there are still differences between the 
median flux levels between quarters. To mitigate the effects of this 
we normalise the light curve. A representation of the true flux level 
was obtained by first normalising independently for each quarter: 

_ Fi — med (F) 



MAD (F) 



(1) 



where -F llorm is the vector of normalised fluxes, F is the vector of 
non-normalised fluxes, Fi is the i th element of vector F, med (F) 
is the median of F, and MAD (F) is the median absolute deviation 
from the median (MAD) of F. MAD is defined as 



MAD (F) = med - med (F) 



(2) 



We show the full custom extracted and cotrended time series 
containing seven quarters of data in Fig. [3] The light curve is that 
of a SU UMa type DN (cf. t he adapted lig ht curve of VW Hyi 
from Bateson (1977) shown in Warner 2003). We observe five out- 
bursts (labelled OB 1-5 in Fig. [3}: four NOs and one SO. One of 
the NOs (OBI) immediately precedes the SO (OB2) while another 
NO (OB3) is of much lower intensity than the other NOs, and oc- 
curs a mere 4 days after the end of OB2. Our method of normali- 
sation is relatively crude; if, for example, there are significantly 
different levels of contamination or the noise level was dramat- 
ically different between the quarters, we would be incorrectly 
estimating the relative outburst amplitudes. We gain confidence 
that the normalisation is appropriate for the following reasons; the 
outbursts in Q5 an Q7 are of comparable intensities, and there is a 
relatively smooth transition from the tail of the SO seen in Q3 and 
into Q4. It is possible that we are severely underestimating the am- 
plitude in Q4 and OB3 is actually of comparable intensity to OBI, 
OB4 and OB5, but this would require a sharp break in the decline 
of the SO in the data gap between Q3 and Q4. 

2.4 Archival Observations 

We were unable to locate NIK 1 in quiescence in Kepler full-frame 
images which have a limiting magnitude of Kp ~ 20, and so we 
looked to archival data to find NIK 1 in a quiescent state. Palo- 
mar Sky Survey I and II ( POSS-I and II) obs ervations appear in 
the USNO-B1.0 catalogue jMonet et al.ll2003l) . USNO-B1.0 lists 
a source offset (in the direction of KIC 4378554) from NIK 1 by 
2.7 arcsec. The POSS-I blue magnitude of this source is given as 
19.15 and the source is catalogued as USNO-B1 1294-0347652. 
They again find this source in the POSS-II blue image and derive a 
blue magnitude of 15.13. 

POSS-I and II images are available of this region as part of 
the Digitized Sky Survey (DSS). We examined the POSS-II images 
of this region and found a source we estimate to be at B ~ 17.3 
based on the brightness of nearby stars in the blue image, much 
brighter than the detection limit of the full-frame images where we 
do not detect NIK 1. We deduce that this source is NIK 1 in a state 



of outburst. Fig. [4] shows a comparison between the POSS-II red 
and blue images which were taken at different epochs. No source 
is detected at the position of NIK 1 in the red image. 

Despite the USNO-B1.0 listing a source at the same position 
in POSS-I blue data we were unable to locate any source. However, 
the fact that the source is listed with two quite different bright- 
ness measurements leaves us in little doubt that USNO-B1 1294- 
0347652 and NIK 1 are one and the same. 



3 OUTBURSTS 

The Kepler data contain observations of five outbursts: four NOs 
and one SO. The dates the outbursts reach maximum and their rise 
and fall durations are given in Tablef3] 



3.1 Normal outbursts 

The four NOs are shown in detail in Fig. [5] and [6] OB4 and OB5 
are very similar looking with comparable brightnesses. However, 
we caution that the uncertainty in measuring the outburst peak flux 
is high because of the method used to normalise these data. The 
sha pes of the NOs a re typical of those seen in other DN (e.g. V344 
Lyr, ICannizzo etai]|2010l) consisting of a fast rise to peak bright- 
ness of between 0.81-1.25 days with a somewhat slower fall to 
quiescence taking 1.86-2.57 days. 



3.2 The superoutburst 

The SO we detected, OB2, is shown in Fig. [5] It reached peak flux 
on BJD 2455176.86 and returned to quiescence 10.42 days later. 
We note that the SO was immediately preceded by a NO, OBI, 
which did not return to quiescent brightness before the onset of the 
SO. OB3 then occurred only 4 days after the SO returned to the 
quiescent flux level. 

3.3 Superhumps 

Superhumps always o ccur during SOs (e.g. IVogtlll980l : rWarned 
l2003l ; lKato et alj|201 ll) . In NIK 1, the superhumps start during the 
decline of the NO, OBI, at the stage when the SO is likely begin- 
ning and continue until the OB3 begins. 

We performed a Fourier transform of the SO searching for pe- 
riods ranging from the Nyquist period (58.84 min) to 12 h. The 
strongest frequency (Fa) found is at 13.03 ± 0.03 cycles per day 
(1.842 ± 0.004 h) whe re the uncertainty wa s calculated using a 
Monte Carlo technique dLenz & B reger 2005). In addition, we de- 
tect a signal at a frequency of 22.85 ±0.08 cycles per day (=1.05 h, 
Fb)- The first overtone of the strongest frequency should occur at 
26.06 cycles per day (Fi). However, because the Nyquist frequency 
(Fn) of the data is 24.47 cycles per day, we cannot uniquely de- 
tect Fi- The value of F 2 - F N is 1.59 ± 0.06 and F N - F B is 
1.62 ± 0.08, therefore, we concluded that the frequency at 22.85 
cycles per day is very likely to be a reflection of the first overtone 
about the Nyquist frequency. 

In order to better characterise the superhumps, we propagated 
the time series through a boxcar high-pass filter with a cut-off fre- 
quency of 18.0 d _1 . This is shown in the top panel of Fig. [7] We 
calculated a discreet Short Time Fourier Transform (Harris| [l978l) 
over the duration of the SO. We used the filtered data and calculated 
a Fourier transform every 0. 1 d with a siding window size of 2.0 d. 
The resulting power spectrum is shown in the lower panel of Fig|7] 
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Figure 3. The fully processed light curve of NIK 1 showing seven quarters of data. The flux time series has been extracted from the target pixel files, corrected 
by fitting and subtracting the CBV and normalised with weights. The normalisation was done by first dividing each quarter of data by its median and then 
dividing this value by the median absolute deviation of the corrected time series. We then divided the whole data-set by the median flux of that data. The five 
outbursts we observed are marked with OB 1-5 and beneath the time series is shown the quarter over which the data were taken. 
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Figure 4. The DSS2 blue image taken on JD 2447351.8 is shown in the left panel and the DSS2 red image taken on JD 2448508.7 in the right panel. Circled 
in white on the right panel is a source in the same position as we identified NIK 1 in the Kepler images. A white circle in drawn in the same position in the 
left panel but no source is detected. The green markers are the positions of sources which are in the Kepler Input Catalogue, no source exists in the catalogue 
at the position where we have detected the outbursting star. 
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Table 3. Outburst occurrences. 



Outburst Quarter Type of outburst BJD of F max T r i se Tf a u 



1 


3 


Normal 


2455173.85 


1.08 


1.96 (a) 


2 


3 


Super 


2455176.86 


1.05 (b) 


10.42 


3 


4 


Normal 


2455190.10 


1.25 


1.86 


4 


5 


Normal 


2455342.84 


0.86 


2.57 


5 


7 


Normal 


2455523.17 


0.81 


2.39 



(a) The fall time is truncated by the onset of the SO 

(b) The rise time is truncated by the tail of the first NO 




BJD - 2454900 

Figure 5. The corrected and normalised flux time series of NIK 1, in the region around OBI, OB2, and OB3. The SO is preceded by a NO, and 3 d after 
the SO we see a weaker, NO. There is a gap in the data covering the SO from BJD 2455182.5-2455185.4 due to the quarterly roll of the spacecraft and data 
downlink. 



Both Fa and Fb are visible for the entire duration of the SO and 
persist through the interval between OB2 and OB3. The superhump 
falls below a 3<r detection level at BJD 2 455 190, close to the time 
OB3 begins. We do not detect the two signals over the rest of the 
time series with any significance. 

In the upper two panels of Fig. [8] we zoom in around the two 
strongest frequencies seen in Fig. |7J The peak frequency of both 
the fundamental and reflection of the first overtone frequency vary 
over time. The crosses show the position of the maximum power 
in the sliding window power spectrum. In Q3, the period of the 
main peak and the reflection of the frequency of the first overtone 
appear to move in opposite directions. This is expected if the higher 
frequency peak is indeed a reflection about the Nyquist frequency. 
In the lower panel of Fig.[8]we show how the amplitude of both the 
fundamental (red crosses) and the reflection of the first overtone 
(blue crosses) change over time. The second harmonic frequency 
decays slower than the fundamental. The fundamental frequency 



falls below the 3-cr detection level 0.6 d before the first overtone 
frequency. 

We see changes in the superhump waveform when compar- 
ing the modulation at mid-outburst (top panel in Fig. [9), at the tail 
end of the SO (middle panel), and when the system in a fainter 
state before the rise of OB 3 (lower panel) . The third time series 
is double-peaked which explains the reflection of the second har- 
monic decaying slower than the fundamental frequency in Fig. [8]- 
the waveform now has twice the frequency of the initial waveform. 



4 DISCUSSIONS 

4.1 The nature of the outbursts 

NIK 1 i s now the fourth CV which has been seen in outburst b y 
Kepler dStill et alj[20icl ; iKato et al.|[2oTTI : iGarnavich etai]|201 ll) . 
All four have undergone a SO. In two of the three other cases 
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Figure 6. The corrected and normalised flux time series of NIK 1, zoomed in on the region around OB4 in the left plot and OB5 in the right plot. The large 
data gap in the OB5 plot is owing to the monthly Earth-point occurring. No data is recorded during this time. 



the SO has been immediately preceded by a NO dStill et alJ|2O10| ; 
iKato et al.l201 1). In the thi rd case, the rise of the S O is missed dur- 
ing a telescope safe-mode dGarnavich et alj|201 ll) . While statistics 
are still low this raises a question: are all SOs triggered by a NO? 

iBatesonl d 1977b sugges ted that in the DN VW H yi, some SOs 
are triggered by a NO, and iMarino & Walked JT979) posited that 
this is the case for all SOs of VW Hyi. In a few other ground-based 
observations of D N, the SO ha s been obse rved to be precede d by 
a NO - e. g. T Leo dKatolll997h. TV Crv dUemura et alj|2005h and 
GO Com llmada et alj2004l) . ICannizzo et alJd2o7(il) postulates that 
a NO occurs due to thermal instability in the disc, and, when there 
is enough mass built-up from previous outbursts at large radial dis- 
tances, a SO can be triggered. However, it takes time to propagate 
the heating front to large radii where this over-density of material 
in the outer disc exists. We first see the NO and then when outer 
material transitions to the hot state, the SO occurs. 

The limit cycle disc instability model predicts the disc is much 
smaller directly after a SO and easier to excite into an outburst state. 
The intensity of outbursts should be lo wer and outbu rsts should oc- 
cur more frequently for a smaller disc dLasotal200lh . Our observa- 
tion are consistent with this pattern. We observe no outbursts before 
OB 1 - the system was quiescent for 208 days - and in the remain- 
ing observations the longest gap between outbursts is 181 days, the 
gap between OB4 and OB5, the final two outbursts. The NOs are 
more frequent after the SO and the spacing increases over time. 
We also see a much smaller outburst directly after the SO, which is 
consistent with theory. 

We observed superhumps lasting for ~ 17 d. They begin 
during the fading of the NO which preceded the SO and last un- 
til the beginning of OB3. We see superhumps for 3^t days af- 



ter the SO has faded to quiescence. There have only been four 
othe r reports of superhump persisting a fter the end of a SO, VI 159 



other reports ot supernump 
Ori dPattersonetal.lll 995). 



„, ER UMa ( Gaoetal. 
dPattersonetalJll998h . and V344 Lyr dStill et al 



199' 



201 



EG Cnc 
It seems 

likely material will persist at large radii well after the majority of 
material has been accreted by the white dwarf and the reason su- 
perhumps have been seen persisting so few times is because of a 
lack of sensitivity in ground-based observations. 



The disc instability model dOsakill 1 9891 : lLasot J200 lh predicts 
two main sources driving the superhump modulation; viscous dis- 
sipation within a flexing disc driven by the compression of the 
disc opposite the secondary star are single peaked, and the hot- 
spot where the accretion stream im pacts the disc swee ping around 
the rim of a non-axisymmetric disc dWood et"ai]|201 ll) . During the 
early superhumps, the hot-spot signal is much weaker than the sig- 
nal from the disc and the phase curve is single peaked. As the disc 
fades, the hot-spot from the impact of the accretion stream onto the 
disc will become brighter relative to the disc because the disc has 
been drained of material and so is much less luminous. 

We see the sumperhump waveform change over the course of 
the SO. In the folded time series shown in the upper panel of Fig. [9] 
the waveform is asymmetric, the pulse rises around twic e as fast as 
it fall s. The shape is very similar to that seen in V344Lvr dstilletal.l 
l2010h at a comparable t ime in the SO evolu tion and in the simula- 
tions shown in fig. 1 of IWood et all d201lh when modelling early 
superhumps. In the middle panel, we see some evolution of the 
waveform and no longer see evidence f or the asymmetry; t he light 
curve is similar to that seen in fig. 3 of IWood et al] d201 lh where 
they have accretion stream impact in the model. Our assumption is 
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Figure 7. The upper panel shows the SO after being propagated through a high-pass filter. The lower panel shows a Short Time Fourier Transform of the 
filtered data shown in the upper panel. Two periodic signals are detected: one at 1.842 h, the other at 1.05 h. These signals persist into the Q4 data and can be 
seen after the NIK 1 has faded before OB 3 occurs. 



that in the middle panel we are seeing no variable component from 
the accretion stream. 

The waveform changes in the lower panel of Fig. [9] Similarly 
to the other panels there is a peak around phase 0.9 - any difference 
can be attributed to the error on the period on which the data we 
folded. However, there is also a secondary peak which is offset in 
phase by ~ 180° and of lower inte nsity than the main peak. This is 
wholly consistent with the fig. 2 of lWood et al.l d201 if) where both 
the disc and the hot-spot contribute to the emission for the system. 

Given that orbital peri ods are typica lly a few percent lower 
than the superhump period JWamedl2003T) . we can place an upper 
limit on the binary orbital period of 1.842 ± 0.004 h, this places 
NIK 1 below the period gap. This meets expectations as a lmost all 
SU UMa stars are below the period gap dKato et alj201lh . 



4.2 DN in the Kepler field of view 

What is the possibility of finding more DNe as ba ckground source s 
in the Kepler data? Using the empirical relation o f |Pattersodj2rJTTh . 
the lowest peak absolute magnitude of a dwarf novae outburst in a 
system with a period longer than 90 min is My = 5.3. This im- 
plies Kepler can observe DNe in outburst out to a distance of 8.7 
kpc if we assume we can observe all sources with a peak bright- 
ness Kp < 20, which is approximately the source confusion limit 
of Kepler. However, we cannot assume an isotropic distribution out 
to this distance. Instead, we restrict our sear ch to only source s less 
than 300 pc above the Galactic plane, which Patt erson! ( 1201 II) esti- 
mates to be the scale height for DNe. 

The centre of the Kepler field of view is 13.5 deg above 
the Galactic plan e and the field of view is 116 square degrees 
dKoch et alJlioTch . Therefore the greatest distance an observable 
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Figure 8. The discreet Short Time Fourier Transform showing the two peaks identified in Fig.[7]in greater detail in the upper two panels. The red crosses show 
the position of the peak in the power spectrum in the region shown. The two panels' colors are scaled the same. The lower panel plots the peak power for both 
the fundamental frequency (red) and first overtone (blue). The black line is the 3-<r detection threshold. 



dwarf nova can be from Earth and still be within our scale height 
constraint is 1285 pc, and our search volume is 2.5 x 10 7 pc 3 (this 
includes the entire field of view, not just the pix ies downloaded). 
The period luminosity relation of Patt ersonl ( l201 ll) predicts all DNe 
outbursts in our search volume will be brighter than V = 15.8 and 
easily be detected by Kepler, even if they have high reddening. 



IPattersonl d201lh gives a lower limit of 1.8 x 10 6 DN pc -3 
but this value is calculated using only 42 systems with known 
distance and most likely under-predicts the true value signifi- 
cantly. Other es timates from observa tions include 7.9 x 10 ~ 6 pc~ 3 
dAketalJl2008l). 1.1 x 10~ 5 pc~ 3 iPretorius et al]|2007t). 2.1 x 
10~ 5 pc~ 3 Mertz et al.lll990l) and 10~ 4 pc~ 3 iShara et alj[l993h . 



Theoretically the number is thoug ht to be around 1.8 x 10 4 pc 3 
dKolbl f 19931 : iPretorius e t al. 2007). For each of these estimates we 
have calculated the inferred number of DNe potentially visible to 



Kepler. These are given in Tableland assume an isotropic distri- 
buti on within t he sea rch range of one Galactic disc scale height. 

Ak et al. (2008) measures the scale height of the dwarf nova 
population to be only 150± 18 pc. If we use this value as the ceiling 
of our search the numbers predicted decrease by a factor of 8. These 
values are also shown in Table[4] 

The predicted values range from 6-4485 which leads us to 
believe that the true space density of dwarf novae is wholly uncon- 
strained. 

If we are considering systems that may be found as back- 
ground sources in apertures of target stars, the numbers decrease 
by a factor of ~ 17 because Kepler only downloads the pixels from 
5.75 per cent of the focal plane. With this restriction, the predicted 
number in Kepler ranges from 0.3-258. With only 9 known DNe 
in the data we are unable to rule out any of the models at the 3-cr 
level. The predictions are likely to be an underestimate of the true 
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Figure 9. Three sections of the light curve of NIK 1 folded on a period of 1 .842 h and binned. The top panel covers the end of the last 2.5 days of the Q3 data, 
the middle panel the tail of the SO, which is the first two days of Q4 and the lower panel covers the quiescent period after the SO and rise of OB3. The flux 
levels are arbitrary but are internally consistent between the three plots. 



Table 4. The number of dwarf novae predicted to be within one scale height of the Galaxy and in the Kepler field of view are given in this table. We use two 
estimates of the scale height; 300 pc and 150 pc. Likewise we use six different estimates of the space density which range from 1.8 X 10~ e to 10 — 4 pc — 3 . 





Sources in focal plane 


Sources in downloaded data 


Space Density 


Scale Height [pc] 


Scale Height [pc] 


DN pc" 3 


300(1) 


150 (2) 


300(1) 


150 (2) 


1.8 x 10~ 6 (1) 


44.9 


5.6 


2.6 


0.32 


7.9 x 10- 6 (2) 


196.8 


24.6 


11.3 


1.4 


1.1 x 10~ 5 (3) 


274.1 


34.3 


15.8 


2.0 


2.1 x 10" 5 (4) 


523.4 


65.4 


30.1 


3.8 


10- 4 (5) 


2492 


311.5 


143.3 


17.9 


1.8 x 10~ 4 (3,6) 


4485 


560.7 


257.9 


32.2 



(1) Patterson (2011) . (2) AketaL (2008) . (3) Preterms et al. (2007) 
(4) Hertz et al. (1990) , (5) Shara et al. (1993) , (6) Kolb (1993) 
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number because this estimate would not include the source in this 
work as it would be further than 300 pc from the Galactic equator. 
A systematic search of the target pixel files may lead to a bounty of 
new DNe. Such work will help constrain the true value of the space 
density of DNe. 



5 SUMMARY 

We have serendipitously observed outbursts of a DN which was 
captured within the aperture of a G-type star. We first extracted the 
flux of this source from the target pixel files and then removed sys- 
tematic trends from the data using cotrending basis vectors. We ob- 
served four NO and one SO over 22 months. The SO was likely trig- 
gered by a NO. We observed superhumps during the SO and saw 
the waveform change from being single peaked to double peaked 
over the course of the outburst which we attribute to the appearance 
of the accretion stream impact onto the disc. 

Given this DN was observed as a background source we used 
various predictions of DN space density to predict the number of 
DN in the Kepler field of view and in data already collected. The 
models vary wildly from predicting 4485 and 257.9 in the field of 
view and in collected data, respectively, to 2.6 and 0.32. A detailed 
study of the background pixels collected by Kepler could be used 
to constrain the space density of DNe. 
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